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Abstract 

We construct a class of spin foam models describing matter coupled to grav- 
ity, such that the gravitational sector is described by the unitary irreducible 
representations of the appropriate symmetry group, while the matter sector is 
described by the finite-dimensional irreducible representations of that group. 
The corresponding spin foam amplitudes in the four-dimensional gravity case 
are expressed in terms of the spin network amplitudes for pentagrams with 
additional external and internal matter edges. We also give a quantum field 
theory formulation of the model, where the matter degrees of freedom are 
described by spin network fields carrying the indices from the appropriate 
group representation. In the non-topological Lorentzian gravity case, we ar- 
gue that the matter representations should be appropriate SO (3) or SO (2) 
representations contained in a given Lorentz matter representation, depend- 
ing on whether one wants to describe a massive or a massless matter field. 
The corresponding spin network amplitudes are given as multiple integrals 
of propagators which are matrix spherical functions. 
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1 Introduction 



The Barrett-Crane spin foam state sum models of four-dimensional quantum 
gravity [0, ||, represent a promising approach for constructing a quantum 
theory of gravity. These are non-topological state sum models, and they give 
discrete spacetime three-geometry to three-geometry transition amplitudes 

0. One can consider these models as examples of regularized path integral 
approach to quantum gravity. The spacetime is described by a simplical 
complex, and the state sum is over colored triangulations. This approach 
is a category theory generalization of the Regge calculus, where the basic 
category involved is the category of representations of the appropriate sym- 
metry group. This group is taken to be SO{d) in the d-dimensional Euclidian 
gravity case, or SO{d —1,1) in the Lorentzian case. Other groups, like anti 
de-Sitter, can be considered as well 0. 

One can also think of the spin foam models as the higher-dimensional 
generalizations of the string theory matrix models, since the spin foam am- 
plitudes are given by a field theory Feynman diagrams which are dual to 
spacetime triangulations 0, |^, |T^, |Tl|, |^ . In the d = 2 case one recovers the 
string theory matrix models |l^ . 

The standard way of making the spin foam partition function and the 
amplitudes finite is by passing to the corresponding quantum group category 

1, |,|]. However, it was discovered in ||^, ^ that even in the Lie group case, 
the BC model could be made finite by changing the field theory propagator. 



It has been proven that the BC Euclidian model is perturbatively finite |l3 
and recently this has been proven for the BC Lorentzian model |T^, i.e. 
the amplitude for any non-degenerate finite triangulation of the spacetime is 
finite. 

In order to make a spin foam model a serious candidate for a quantum 
theory of gravity, like string theory, one should be able to incorporate matter. 
Mathematically this would mean enlarging the category of representations, 
by adding a new set of representations corresponding to the matter fields. 
The first proposal of this kind was made by Crane [0, who suggested a tensor 
category of unitary representations for the quantum Lorentz group. In this 
proposal a subset of the irreducible representations (irreps) corresponds to 
the gravitational degrees of freedom (DOF), while the rest corresponds to 
the matter DOF. The corresponding state sum model is topological, which 
is justified as an acceptable feature of the theory at high energies. Then one 
conjectures that the local dynamics would appear at low energies, by some 
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sort of symmetry breaking mechanism. In a more recent paper Crane 
has proposed to interpret the graphs of the field theory formulation which do 
not correspond to triangulations of manifolds, as triangulations of manifolds 
with conical singularities, and then to interpret these singularities as matter 
DOF. These are interesting and very geometrical proposals; however, it is not 
clear how to relate these matter DOF to the particle field theory description 
in terms of fields carrying pseudo-unitary finite-dimensional Lorentz irreps. 
In this paper we give a formulation of a spin foam model with matter 



inspired by the particle field theory formalism [0. We consider a class of 



spin foam models describing matter coupled to gravity, where the gravita- 
tional sector is described by the unitary irreps of the appropriate symmetry 
group, while the matter sector is described by the finite dimensional irreps 
of that group. In the non-topological four-dimensional Lorentzian gravity 
case, the gravitational sector will be described by the BC model, while in 
the matter sector we will argue that one needs to take appropriate 5*0(3) 
or 5*0(2) representations contained in a given finite-dimensional Lorentz ir- 
rep. The construction we present here naturally follows from the field theory 
formulation of the spin foam amplitudes f^. It was shown there that a 
spin foam transition amplitude can be represented as an matrix element of 
an evolution operator, where the initial and the final states correspond to 
the three-dimensional spin networks (or spin nets for short) induced by the 
boundary triangulations. One can represent these spin net states as the spin 
net operators acting on an appropriate Fock space vacuum. One can think 
of these spin net operators as the second quantized spin net wave-functions 



from the canonical loop gravity [|18j. Since the spin net wave-functions car- 
rying fermionic matter have been introduced in the context of canonical loop 
gravity it is not difficult to generalize these results to the second 

quantization formalismjj], so that the matter insertions at the spin net ver- 
tices become appropriate creation/annihilation operators. 

In section 2 we give basic definitions and formulas concerning the spin 
net functions. In section 3 we give basics on spin foams, and explain how 
the spin foam amplitudes are related to the spin net amplitudes. In section 



^The spin network quantum field theory is an example of what is traditionally called 
a third quantization of gravity, because it allows for more then one universe (a boundary 
triangulation) in the formalism . Since general relativity is a classical field theory, then 
the expression "second quantization" is a more appropriate terminology. However, the 
exact relationship between the loop gravity formulation and the spin-foam formulation is 
not yet fully understood. 
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4 we construct the fermionic spin net operators and states, and describe the 
general structure of the interaction with the gravitational background. In sec- 
tion 5 we discuss the fermions in the Euclidian gravity case. We show that 
in the non-topological gravity case, the corresponding spin-net amplitudes 
can be written as the multiple integrals of propagators over the homogenous 
space = S0{4:)/ S0{3), where the fermionic propagator is given by a ma- 
trix zonal spherical function. These results are straightforwardly generalized 
to the Lorentzian gravity case in the next section. In section 7 we discuss 
the general structure of the interaction terms for the case of matter fields 
of arbitrary spins. In the non-topological gravity case, the spin-net matter 
propagators are described by the corresponding matrix spherical functions. 
We also argue that the spin-net propagators for the analogues of the mass- 
less gauge bosons are the matrix spherical functions restricted to an 5*0 (2) 
subgroup. In section 8 we present our conclusions. 

2 Spin net functions and amplitudes 

We review the results on spin nets from a point of view which is most 
suitable for the type of generalizations we want to make. 

Consider square-integrable functions on G^, where G is a Lie group 
and E is a. natural number. A gauge invariant function is determined by 
an oriented graph F with V vertices and E edges, where a group element is 
associated to each edge e & E. If we label the vertices of T as Vj, 1 < j < V, 
then the edges connecting the vertices j and k can be denoted as ejk, and 
the corresponding group elements as gjk- Orientation is such that gjk = g^l- 
A function \E'r associated to the graph F is gauge invariant if 

^r{9jk) = '^r{hjg,kK^) , (1) 

for all hj,hk G G. 

By using the Peter- Weyl theorem, one can show that a gauge invariant 
function \E'r can be expanded as 

(5-1 , • • • , S'b) = CA,i ^^,A,^ i.9u---^9E) (2) 

where A = {Ai, A^;} denote the unitary irreps of G labeling the edges of 
F, and i = {ii, iy} denote the intertwiners associated to the vertices of F. 
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The basis functions $ are called the spin net functions, and they are given 

by 

U D^Klige) U c::}''^ . (3) 

eS-B veV 

where D^^\g) is the representation matrix of gf, acting in the space V{A). All 
the representation indices cce and (3e are contracted by the indices from the 
product of the intertwiner components. Contraction is defined as X'^Ya = 
C^^X^a. where C"'^ are the components of the V{k) V{h*) = V*{A) 
intertwiner, and V* is the dual vector space. Since any intertwiner can 
be represented as a linear map V^(Ai) ® • • • ® V(A„) — >• C, we denote the 
corresponding vector space as Inv{Ai, A„). 

Given a spin net function $r,A,i(5'e), one can obtain a group invariant 
associated to the spin net 

A{r,A,i)^^r,A,ii9i^9o,--;9E^go) , (4) 

where go is the identity element of G. We will call the number A the am- 
plitude for the spin net (F, A, i). These objects are the basic building blocks 
of the state sum amphtudes. Since D^{gQ) — S^, the amphtude A is simply 
the product of the intertwiners C'"*^*"-' for all vertices v E V, with all the 
representation indices being contracted. The graphs F which are relevant for 
the state sum models are the one skeletons of (i-simplices, as well as the 0,^ 
graphs (two points joined by d edges). 

In the case of open spin nets, the above results can be easily generalized. 
The amplitude A for an open spin net will be a group tensor given by the 
product of intertwiners for the vertices, where all the representations indices 
are contracted, except the indices corresponding to the irreps labelling the 
external edges. This can be represented as 

A^^];;- (F, A, i) = C;- ^''^ ■ ■ ■ C-—;:^^ ■ ■ ■ C;;f ^^'^ • • • C;;; ^'^^ . (5) 

Hence the amplitude for an open spin net is not an invariant, but it is a 
tensorial quantity. One can construct a group invariant by contracting the 
external edge indices with the vectors and the co-vectors from the external 
irrep vector spaces 

S{T,A,i,v,u)=Ai\::\T,A,i)v0,---u"^--- , (6) 

where are the components of a vector v € ^(A), while -u" are the compo- 
nents of a covector u e V*[A). This formula will be useful for constructing 
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the invariant actions, since it is the group theory description of the actions 
used in particle field theories. 

It is now also clear how to construct the open spin net functions (ge)- 
Note that the expansions of the functions over G"^ used in the field theory re- 
formulations of the d-dimensional spin foam amplitudes as Feynman graphs 



@, ^, [T^, |TT], |], are simply expansions in terms of open spin net functions 
associated to the graphs {d edges emanating from a single vertex). 

In the case of the BC model [|l], 0, the amplitude for the x graph is not 
just the intertwiner C^^ G Jnf( Ai, A4), but it is a linear combination of 
these 

5x=E«^^x^ • (7) 

i 

This is the consequence of the fact that the subcategory of irreps which are 
used have invariant vectors under the 5*0(3) subgroup. 



3 Spin foam amplitudes 

The state sum models which have been used for constructing (i-dimensional 
quantum gravity models, for a review see ^ , are based on a two-complex 
J = {V,E,F), consisting of a finite number of vertices v E V, edges e E E 
and faces f E F. The two complex J is taken to be dual to a triangulation 
of the spacetime manifold M. Since the faces of J are labeled with the irreps 
of G, one can regard this object as a generalization of the spin net, where 
the one-complex F = (V, E) is replaced by a set {E, F) from the two complex 
J. That is why J is called a spin foam pSf. Hence these state sum models 



give the amplitudes for the spin foams. An amplitude for a spin foam can be 
written as 



A /eF e&E veV 

where A2 is the amplitude for a face, Ai is the amplitude for an edge, and 
Aq is the amplitude for a vertex of the two complex J. 

There are several ways to determine the amplitudes An, but what is 
important for our purposes is the fact that these amplitudes are proportional 
to the amplitudes for the closed spin nets defined in the previous section. In 
the d = 4 spacetime case, Aq is the amplitude for the 4-simplex graph (a 
pentagram) |l|, |^. Ai can be one, or it is the amplitude for the 84 graph 
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in the modified BC model P, 0, IH. A2 is usually taken to be diniA, which 



is a single loop spin net amplitude. One can also take A2 to be any function 



of A, such that a group Fourier transform of A2 exists pO[| . 

Another aspect of the spin foam amplitude (H) which will be important 
is the fact that when A2 = dimA then the amplitude can be represented as 
a Feynman diagram of a field theory over the group G |, 0, |ll|, ^. The 
action for this field theory can be written in terms of the Fourier modes (j), 
which transform like tensors from Inv{Ai, ...,Ad) spaces. In the 4 = 4 case 
the action takes the form 

A,i 

+ ^E^"i--5(xi,---,X5)CV--r4 , (9) 
where 0" = 0A^...A'^(j)- The vertex functions Q and V are given by 

Q-l(Xi,X2) = 5Mi(Ai,---,A4;^l,22) (10) 

V(xi,---,X5) = 5^%(Ai,---,Aio;zi, ...,^5) , (H) 

where 5 are the appropriate delta functions of the indices. Ai can be the 
amplitude for the 64 graph, or it could be one, when is the amplitude for 
four parallel lines. Aq is the amplitude for the pentagram. In the case of BC 
models, Aq and Ai will not depend on the intertwiner labels i^- 

The corresponding Feynman diagram is just the one complex (V", E) ob- 
tained from the two complex J = (V, E, F). This diagram is a {d+ l)-valent 
graph in the case of d-dimensional model. Its value is given by (H), which can 
be represented by a collection of spin net diagrams, such that to each vertex 
we associate a cZ-simplex graph, while to each edge we associate d parallel 
lines, or a 0^ graph in the modified version. 



4 Fermionic matter 

The Fourier modes 0" can be promoted into creation and annihilation oper- 
ators 0, so that the spin foam amplitude can be represented as an matrix 
element of an evolution operator, in analogy with the particle field theory 
case. As shown in [|] , one can construct the "in" and the "out" states which 
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describe the spatial spin nets corresponding to the boundary triangulations 
of the spacetime manifold. Such spin net states can be represented by 



|7)= n 4 JO) , (12) 

dGV(7) 

where |0) is the vacuum of the Fock space constructed from the creation and 
annihilation operators and 0^, and 7 is a four-valent spin net graph dual 
to the boundary triangulation. All the representation indices in ( [T^ ) are 
appropriately contracted. 



In the loop gravity formalism ||21|, the fermionic matter is introduced 



by replacing the spin net 7 with an open spin net 7ir, which is obtained 
from 7 by putting an external edge carrying a fermionic irrep at each site 
of 7 where a fermion is located. One can also put more than one fermion 
at each lattice site, up to dimAg (Pauli exclusion principle). In that case, 
the external edge is labeled with an irrep from an antisymmetrized tensor 
product of k Ag irreps |p2 |. 



In the spin foam formalism, the analogue of the loop gravity construction 
would be to construct the |7i?) state by introducing the fermionic creation 
and annihilation operators iptri^v) and '^^(x^), where a is the representation 
index of A^, and the label x„ denotes the four irreps of the spin net site where 
the fermion is located. Hence the corresponding state would be given as 

M= n n 4jo), (13) 

where V is the set of vertices where the fermions are located. 

In order to calculate the transition amplitudes between two fermionic 
spin net states, which would be a fermionic spin foam amplitude, we need 
to introduce the corresponding terms in the field theory action (^. Since 
classically the fermionic propagation is described by a line in the spacetime, 
and since the fermions interact with the gravitational field, the candidate 
action will be 

'^F = EV^''(x)^'''(x')V..'(x,x',Xi,---,x„)0x,---0x„ + (/i.c.) , (14) 

A 

where {h.c.) stands for the hermitian conjugate term. 

The form of the vertex V can be restricted by two requirements. First, Sp 
has to be a group invariant, i.e. a scalar. Second, we require that the fermions 
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propagate only on complexes which are dual to the Feynman graphs of the 
gravitational model @. The idea behind this is that the purely gravitational 
model builds up the space-time on which the fermions propagate. Hence 
n = 5, and 

Vaa'{x, X', Xi, ■ ■ -, X5) = (5^°(5x,Xi^x',Xfe^aa'(Al, Alo) , (15) 

where Acra' is a spin net amplitude for the pentagram with two external edges 
attached at site 1 and site k, where 1 < A; < 5. 

Note that one can take a more complicated spin net graph for the ampli- 
tude Ao-o-', by adding internal edges labelled by the matter irrep, see sect. 5. 
One would then need to study the semiclassical limit in order to determine 
which spin net would give a correct physical theory. Since that is beyond the 
scope of this paper, and for the sake of simplicity, we will analyze only the 
simplest possible spin nets. 

We will also add to Sp a purely quadratic fermion term 

j2r{x)Q..'{x,x')r'ix') + ih.c.) , (16) 

A 

in order to have a well-defined perturbative expansion. The propagator 
will be determined by the amplitude for the 64 graph with two external edges 
labelled with A^, or it could be a trivial amplitude for five parallel lines. We 
denote the propagator amplitude as Gaa'- 

Note that (ipa)* = ip'^ = C"" ip^', and the vector space duality * is in gen- 
eral different from the complex conjugation (reality) properties. The reality 
properties determine relation between the creation and annihilation opera- 
tors ip and ip'^ [^]. 

Therefore the maximal closed spin net subgraph for Acra' has to be at 
least the pentagram graph. Hence a relevant fermionic Feynman diagram Tp 
which is generated by the 5*4 + Sp action will be determined by a spacetime 
skeleton diagram F (a 5-valent graph dual to a spacetime triangulation) plus 
a line connecting the vertices of F. Therefore Tp will be a graph consisting 
of 5-valent and 7-valent vertices, see Fig. 1. Out of these diagrams we will 
consider only those for which the fermion path is a line which connects the 
centers of a string of adjacent 4-simplices, extending from the initial to the 
final boundary. 

Hence the fermionic spin foam Jp will be given by the usual spin foam 
J and a line L of edges starting from an "in" edge and terminating with an 
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"out" edge. The corresponding amplitude will be given by 

= E n MmiiiL) n ^e(A) n m^) , (17) 

A feF e£E' v£V' 

where 

A^dL) = A:1{v^) G^lieu) A^iv^) ■ ■ ■ G^l^e..^,,) A^Hv,) . (18) 

The number of the "in" and the "out" fermions has to be the same for a 
free theory, and the above formula can be then easily generalized to the case 
when there are several lines L. 

The amplitude A for a fermionic vertex of Jp will be given by a spin net 
amplitude for a pentagram with two external fermionic edges, Fig. 2 

A^,(Ai, ...,Aio)=EC<x'xC^CxCxCx , (19) 

a 

where Cq-'xjC'x ^ /nf (A^, Ai, A4) and Cx G Jnf (Ai, A4). The ampli- 
tude for an fermionic edge of Jp will be given by the propagator amplitude 
, which in the non-trivial case is given by the amplitude for the G4 graph 
with two external edges. Fig. 3 

G^,(Ai, ...,A4)=EC.'xC^x • (20) 

a 

In some of the cases the minimal amplitudes (|T^) and ( PP] ) are zero, be- 
cause Cax = 0. Then one can replace it with C^a'x from Inv{As, As, Ai, ... , A4), 
so that 

A:,iA„...,A,o)=J2C:'xC^xC^C^C, . (21) 

a 

The corresponding spin net graph is a pentagram with two external and 
one internal fermion edge, joined in a line. Fig. 4. Similarly, for the edge 
amplitude , we will have in the non-trivial case the 65 graph with two 
external edges. Fig. 5, so that 

G^,(Ai, ...,A4) = EC<xVxC7 . (22) 
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5 Euclidian case 

We now study the four-dimensional Euclidian gravity spin foam model [||, ^, 
so that G = SO (4). The irreps of SO (4) can be labeled as A = {j,k) 
where j and k are the half-integers labelling the irreps of SU{2). In the 
topological gravity case, one uses all 5*0 (4) irreps for labelling the triangles 
of a triangulation of M. The fermions will be represented by Aj" = (|,0) 
and A~ = (0, |) irreps. These irreps correspond to chiral, or Weyl fermions, 
while a Dirac fermion would be a reducible representation (i, 0) © (0, |). 

The relevant intertwiners for constructing the fermionic spin foam ampli- 
tude are from the space Jnt>(As, Ai, A4). A basis \i) is determined by 



In order to simplify the notation, we take XaY^ = X'^Ya = C^^XcXp- Note 
that the formula (^) can be represented graphically by Fig. 6. Hence 
the intertwiner label is given hj i = (A, A'), and we can now construct 
the fermionic amplitudes for topological gravity spin foam model by using 
the procedure described in sect. 4. However, we will not discuss further 
the topological gravity case, and we will concentrate on the non-topological 
model, because of its importance for physics. 

The set of relevant irreps in the non-topological gravity case is given by 
the simple irreps N = In that case C^^^^'"^^ = 0, which can be proven 

by using a Z2 grading of the 5*0 (4) irreps. This grading is introduced by 
splitting the irreps into even and odd, or bosonic and fermionic irreps via the 
parity function p {n) = { — 1, 1} if n is odd, even, respectively, as 



Hence the simple irreps are even, i.e. bosonic, while A^ are odd, i.e. fermionic. 
Since p (Ai (S) A2) = p (Ai)p (A2) = p(A) where A is an irrep from the tensor 
product, then a product of bosonic irreps can never yield a fermionic irrep, 
and hence (^^^^i- -^^ _ q Therefore the minimal option for the intertwiner 
space is Inv{As, A^, Ni, A^4), and a basis \i) can be constructed from 




(23) 



P U, k)=p {2j)p {2k) 



(24) 



ihs^sNi...N4,(i) _ /^KshsBifyNiN2B2r-iBiB2Bz(^BzNs,Ni 
(7(T'ai...Q4 ~ ^au'Pi ^aia2l32 '-'/3i/32/33 '-"/33a3«4 

/3i,/32,/33 



where B are the bosonic irreps. The equation (p5| ) can be graphically repre- 



sented by Fig. 7. 
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Now we have to impose the constraints coming from the fact that the 
simple irrep vector spaces have invariant vectors under the S0{3) = SU{2) 
subgroup. This amounts to constructing an invariant vector in the corre- 
sponding Inv{n) space. When there are no fermions, then n = 4, and the 
invariant vector, or the BC vertex, is constructed as m 

01. ..04 / y 01. ..04 / y / y ai020 aa3a4 ' \ ^) 

M M a 

where M are the simple irreps. The 3j symbols for 5*0(4) can be expressed 
as the product of 3j symbols for SU{2) as 

/-iNiN2N'i _ ryjlj2j3 f<jlj2j3 (07) 
010203 'mim2m3 nin2n3 " V '/ 

The consequence of using the BC vertex for constructing the spin net 
amplitudes is that these amplitudes can be represented as multiple integrals 
over the coset space X = G/H, which is a three-sphere = S0{4:)/ S0{3) 
25, 26|. For example, the amplitude for the 64 graph will be given by 

e(iVi, iV4) = E sZ\\\-Z'sZ\:f: ■ m 

a 

One can show that 

e{N,,...,N,)= f dxdyllKMXx,y) , (29) 

where 

ir^(x,,)=TrD(^)(x-^.,) = ^^^^M±M , (30) 

sm (7 

and 6 is the geodesic distance between the coset points. This formula works 
because the coset points x and y can be considered as SU{2) group elements 
due to equivalence of and SU{2) spaces. 

We will now use this property to construct the fermionic BC vertex. 
Consider the following integral associated to the 64 graph 

r 4 

Qss'{Ni,...,Ni) = dxdyKss'ix,y)l[KNXx,y) , (31) 
-'^^ i=i 

where 

K,Ax,y)=D^'f\x-'-y) . (32) 
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One can now show that 



where S" is an element of Inv{Ni, A'4, N1/2) given by 

N1...N4N1/2 _ ^NiN2M^'^^^l/2N ^NN3N4 (nA\ 

M,N a,/3 

and A^i/a = 

When ( plD is compared to (|25|) , then it is easy to see that the expression 
( [5^ ) can be also taken as an intertwiner from the space Inv{A'^, A~, Ni, N4). 
due to = A+ ® A~. Therefore the expression (|31|) can be interpreted 
as the amplitude for the open spin net consisting of the 64 graph plus two 
external edges labelled with A+ and one internal edge labeled with A~, Fig. 
5. 

Hence the action Sp for the fermions will be given by 



N 

+ EV^±(x)V^±(x')V<x.(x, x'; Xi ■ ■ ■ X5)0xi ■ ■ -0X5 
N 

+ {h.c.) , (35) 
where and V are determined by the 

G„,{N^...Ni) = f (fxK„,{xi,X2)]\KNAxi,X2) (36) 

and 

^<^^(A^i...A^io) = j ^ d^xK^^{xi,Xk) Y[ KN(ij){xi,Xj) (37) 

spin net amplitudes, respectively. The edge amplitude ( P^D is determined by 
the graph from Fig. 5, and the vertex amplitude (|37|) is determined by the 
graph from Fig. 4. 

The reason why (|31| ) and (|33| ) are the same is that the spinor propagator 
(32) is an example of a matrix spherical function [27|, which is a generaliza- 
tion of the scalar propagator (pOf). In general case if if is a subgroup 
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of G and if r is an irrep of iJ, consider irreps A of G which contain r when 
decomposed with respect to the subgroup if, i.e. 

^A= K'©-=^r©... . (38) 

l<fc<n 

Let ^{g) = PrD^^\g)PT-, where is the projector from Va to if^. ^{g) 
are called r-spherical functions, and P^D^^\g) are eigenfunctions of the 
Laplacian on the homogenous vector bundle Er = {G/H,Vr). Since ^{g) G 
Hom{Hr, Hr), we can construct (j){g) = tr^{g) G Hom{Vr,Vr), where tr is 
a partial trace over the degeneracy label k. Then consider the object 

KA,rix,y) = U{x,xo)(l)ig^^gy)U{xo,y) , (39) 

where f/'s are vector bundle parallel transport operators along the geodesies 
{xo,x) and {xo,y), and gx = <j{x) is a local section of the principal bundle 
{G/H^G). By using g^ = xh, gy = yh, U{x,xo) = D^''\h) and U{xo,y) = 
D^^\h~^), one can show that K is a function of only the coset space points. 
Note that when r is the identity irrep, the expression ( pQ] ) reduces to the 
scalar propagator Kn, where is the class one irrep. Hence (|39|) is the 
matrix generalization of K^, and we can use it to construct the spin net 
amplitudes for graphs with external edges. 

In the 5*0(4) case, let A and B be the generators of the two commuting 
SU{2). Then the 50(4) rotations generators J^,^ can be expressed as 

J = A + B , K = A-B , (40) 

where Ji = ei^^Jjk and Ki = JiQ. The 5*0(4) irreps A are then determined 
by the pairs of half- integers {A, B), where A and B are the angular momenta 
numbers of the two SU{2). We then have 

D^'^'^^e^v) = exp{ieJA + ivKA) 

= exp {ie+AA + id- Bb ) 

= D^^\e+)®D^''\e_), (41) 

where 

e± = e±v . (42) 

The diagonal SU (2) subgroup vectors are given by 

\3,m)=C^^J!\a)®\h) . (43) 
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These vectors form an irrep space of the spatial rotation generators J. In 
the case of A^, we take r = j = 1/2, so that P^- = 1. From (|41|) we get 

K^/,ix,y) = Z}(V2)(^-i . = , (44) 

where 

(^ ^ cos(r/2) J2 + i an sin(r/2) , n = = - , (45) 

and a are the Pauli matrices. 



6 Lorentzian case 

In the Lorentzian case G = 50(3, 1) = SL{2, C), and the finite-dimensional 
irreps (j, k) become pseudo unitary, because the relation becomes 

J = A + B , K = {-i)iA-B) , (46) 

so that the boost generator K becomes an anti-hermitian operator. This is 
the reason why there are no finite-dimensional unitary irreps for the Lorentz 
group. 

Hence the unitary Lorentz irreps are infinite- dimensional, and can be 
labelled as (m, p), where m is a half-intiger, while p is a non-negative real 
number [?]. The gravitational DOF are again described by the simple irreps, 
or the class one irreps with respect to the spatial SU (2) subgroup. These 
are given by the (m, 0) and (0,p) irreps 0. However, only the (0,p) irreps 
are used 0. We could then use the pseudo unitary irreps (j, k) to describe 
the matter fields, but we have to check first the consistency of the state-sum 
model with unitary and pseudo unitary irreps of the Lorentz group. 

We now have = (0, p) in the gravitational sector. The corresponding 
propagator is given by [0] 

' psinh d(x,y) ' 

where d{x, y) is the geodesic distance between the points x and y of the coset 
space X, which is now a hyperboloid 50(1, 3)/S'0(3). The fermions will 
carry the irreps A^. The tensor product decomposition rule of a unitary and 
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a pseudo-unitary irrep can be inferred from the decomposition of these irreps 
under the SU{2) subgroup. For the (m, p) irrep we have [p8[ 



oo 



^(-,P) = . (48) 

j=\m\ 

Since ^(A^) = Vi, we obtain 

H(™,p)®l"(A^)=H(^_i,p)©^(„+|,p) , (49) 
which agrees with |[T^. This imphes 



T^{m,p) ® V{Av) = H(m-l,p) ffi ^(m,p) © ^(m,p) © ^(m+l,p) , (50) 



where A„ = (|, |). 



Hence C^'i'^^^" 7^ 0, and we can now construct the Lorentzian analog of 
the BC spinor vertex (|3^ . For this we need a Lorentzian spinor propaga- 
tor, which is determined by the formula (p9|). It is given by an analytic 
continuation of the Euclidian formula (H), where v —iv, so that 



where now 



K^,2{x,y) = D^^l^\x-^ ■ y) = D^^l^\iv)D^^'^\-iu) , (51) 
D^^I'^\ir) = cos\i{r-l2)l2 + BnsiYi\i{r/2) . (52) 



7 Interactions and gauge fields 

Note that in the discussion in section 4 one could have taken A^ to be an 
arbitrary irrep of G, and moreover, one can take a set of different irreps 
Asi = 5*1,..., Asj. = Sk to represent fields of different spins. In the topological 
gravity case, one can write the action as 

Si = EE^^^(x)^5(x')Q55(x,x') 

S A 
S A 

+ EE^5i(x) ■ ■ ■^5,(x('^))V5i...S,(x, X«; Xl, X5)0xi ■■■0X5 
S A 

+ {h.c.) , (53) 
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where the vertex Q will be determined by the spin net amplitudes for the 
open graphs based on the five parallel lines and the ©4, while the vertex V 
will be determined by the spin net amplitudes for open graphs based on the 
pentagram. A new feature is that the set of possible graphs for the vertex 
'^Si...Sk will contain graphs with more than five vertices, Fig. 8. 

In the non-topological gravity case, we need to implement the H invari- 
ance. This requires a specification of the irreps s of H which are contained 
in the irreps S of G. In that case the matter fields will have the indices from 
the H irreps, and the action can be written as 

S N 

+ J2J2Mx)iJs{x')Vi{x, X', Xl,---, X5)0xi •••0X5 
S N 

+ EEMx) ■ ■■Mx^'W!.:t(x,..., x«; Xl, X5)0 

5 N 

+ (h.c) , 

where the vertices Q and V will be determined by the spin net amplitudes for 
open graphs based on five parallel lines, the ©4 and the pentagram graphs, 
as in the topological case. However, these amphtudes will be now given as 
multiple integrals of the propagators Kj^{x, y) and Ks,s{x, y). Also note that 
the fields ips will now carry the indices from the irreps of the subgroup H. 

In the ETiclidian/Lorentzian gravity case H = SO (3), and it will be in- 
structive to study the case of spin j = 1 matter fields, because in the nature 
the fermions interact through the spin j = 1 gauge bosons. If we take a vec- 
tor irrep 5*1 = (|, |), it will contain a, s = j — 1 irrep of 5*0 (3). We denote 
the corresponding spin net field as Ai{x), while the fermions we denote as 
ijj'^(x) = ipfi^x), where the subscripts ± denote the corresponding 5*0 (4) 
irreps. The spin foam analog of the particle field theory interaction term 
/ d'^xipj^ipA/^ can be written as 

J2^t(x)^|^:'(x')A(x")V'4\x, x', x"; x,, X5)0x, + (h.c.) , 

N 

(55) 

where the vertex V+- is determined by the spin net amplitude for the open 
graph given by the pentagram plus a three- valent matter vertex. Fig. 9. This 
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(54) 



amplitude is given by 



(56) 

Since (|55| ) is a spin foam generalization of the particle gauge field theory 
f/(l) interaction, it is natural to ask what is the analog of the U{1) gauge 
symmetry. One way to deal with the issue of gauge symmetry is to assume 
that the components of the fields appearing in the action ( p^ are the physical, 
or the independent DOF. Note that one can adopt this type of approach in the 
context of fiat space particle field theory [|l7l. One simply takes the creation 
and annihilation operators labelled by the irreducible Poincare representation 
indices and then forms polynomial interaction terms via the formula (|]). The 
relevant Poincare irreps are the massive and the massless ones, so that the 
index labels are the three- momentum, the 5*0(3) irrep indices for the massive 
case, or the two-dimensional Spin{2) = 5*0(2) representation labels in the 
massless case. 

Since one wants the particle interactions to be local in spacetime, one then 
constructs spacetime fields carrying the appropriate Lorentz irreps from the 
Poncare creation and annihilation (c/a) operators, and writes the interactions 
as spacetime integrals over the polynomials of these fields [0]. The problem 
with this approach is that the massless fields cannot be always expressed in 
terms of the helicity c/a operators. This happens for the vector field A^(a;), 
which when expressed in terms of the helicity c/a operators transforms non- 
homogenously under the Lorentz transformations [jT^. However, the field 
f\iu = dif^Ai,] transforms like a Lorentz tensor, so that the cubic interaction 
term would be given by / d^xilj'^^'^ilj f^u- But this is not the type of interac- 
tion realized in nature, which is / d'^xip'-^^ipA^ . This paradox is resolved by 



introducing the U{1) gauge transformations |1T7[ . 

In the case of our spin foam model, this type of problem will not ap- 
pear immediately, because one is dealing only with the c/a operators for 
the Lorentz group and its subgroups. These operators carry the labels 
X = (A''i...A^4) and SU{2) = 50(3) representation indices. There are no 
Poincare irreps labels, and the reason for that is that the particle momen- 
tum is not conserved in curved spacetimes. Hence it is not clear what would 
be a physical interpretation of a spin-foam model based on the Poincare 
group0 

■I- This can be seen on the example of a Poincare spin foam model based on the simple 
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Our assumption that the components of the spin net fields are the phys- 
ical components can be further justified by the fact that these fields live on 
the spin nets defined by the triangulations of three-space boundaries, see 
equation (0). We will then consider these fields as the discretized analogs of 
the reduced phase space canonical formalism fields. The canonical analysis 
of the continuous gravity and matter actions shows that the effect of the 
local Lorentz and the spacetime diffeomorphism constraints is to reduce the 
number of components of the gravitational field, while the matter field com- 
ponents are affected only by the gauge symmetry constraints. The final result 
is that the non-gauge fields have independent components which transform 
as the appropriate representations of SO (3), while the gauge fields indepen- 
dent components transform as the appropriate representations of S0{2). In 
the fiat space limit, these independent component fields become the Fourier 
transforms of the Poincare irreps, so that the non-gauge fields correspond to 
the massive irreps, while the gauge fields correspond to the massless ones0. 

Therefore within this approach the Ai{x) operator cannot be considered 
as a discrete massless gauge field analog, since it carries a j = 1 irrep of 
5*0(3), which corresponds to a massive Poincare irrep. In order to introduce 
a massless-like vector gauge boson in the model, we need a spin net field 
which transforms as a two-dimensional representation of 5*0(2). Therefore 
we need to construct invariant propagators under an 50(2) subgroup of the 
50(3) subgroup. 

Hence consider a subgroup H' of H, and let r' be a representation (could 
be reducible) of H' contained in the irrep r of H. Note that one cannot use 
the K\y propagators for this purpose, where r' is contained in the irrep A. 
The reason is that K\y will be a function on the coset space X' = G/H', 
which is different from the space X = G/H, where the gravitational prop- 
agators K]\f are defined. Therefore these propagators can not be multiplied 
and integrated to give the spin net amplitudes. 

unitary irreps of the Lorentz group. At first sight this appears a natural thing to do, 
because the unitary irreps of the Lorentz group are also unitary irreps of the Poincare 
group, since the Lorentz subgroup is the little group for the zero momentum. Hence the 
gravitational sector would be described by the zero-momentum unitary irreps, while the 
matter sector would be described by the usual particle unitary irreps, having a non-zero 
momentum. However, because the tensor product conserves the momentum, the corre- 
sponding spin-foam amplitude would conserve the momentum of a particle propagating in 
a curved background. 

^Except in the case of massless spin-half fermions, or if there is a Higgs scalar field. 
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Note that the representation r' will induce a vector sub-bundle {G/ H, Vr') 
of the vector bundle {G/H, Vr). We can then define the propagator 

KA,r'ix,y) = Pr'U{x,XQ)(t){g^^gy)U{xQ,y)Pr' (57) 
= Pr,K^^r{x,y)Pr. , (58) 

where Pr' is a projector from Vr to Vr'. It is clear that K G Hom(y' , V) 
and it transforms correctly under x — > xh' and y — > yh' where h' ., h' G H'. 

In our case we take the index on the massless-like spin net field of spin j 
to be from a reducible Spin{2) helicity representation 

V; = C{\j),\-j)} . (59) 

This choice is motivated by the continuum space canonical analysis and with 
the fiat space Poincare invariance. Clearly Vj is a subspace of Vj, and the 
Spin{2) subgroup is generated by the component of the angular momen- 
tum. 

In the j = 1/2 case we have V"(A^) = V1/2 = VY/2) so that the only way 
to distinguish between the massive and a massless fermion is by a choice of 
the terms in the quadratic action (P5|). 

In the j = 1 and Si = (1/2, 1/2) case, the index i in the massless-like Ai 
operator will take two values, so that the cubic interaction term would be 
again given by the formula (0), but now all the indices will belong to the 

5*0(2) Vi and V( representations. 

2 

Note that in the j = I case one can also take = (1, 0) and Si = (0, 1) 
Lorentz irreps. These are the irreducible parts of the antisymmetric tensor 
ffj_iy, and we can denote the corresponding spin net fields as f^{x). The 
spin-foam analog of the action / d'^xip'^'^'^ipffj^i, will be given by 

Y.^f{x)^${x')ft{x")V'^i\x, x', x"; Xi, X5)0x. ■ ■ ■ + {h.c.) , 

N 

(60) 

and the corresponding spin net amplitude is given by 

A%i'= I d''xdyKf{xuy)K\\x,,y)K{'\x2,y)c\;fjX{KM,X^,,xi), 

(61) 

where Ki = i^(i,o),i = ^(0,1), i- 

Note that the construction ( p9D gives for j = 2 a massless-like helicity 
two field h±2{>^)- Since this is a local excitation on a given spin net, it can 
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be considered as a spin foam version of the graviton. Again there are two 
possible choices for the corresponding Lorentz irrep. One can take 5*2 = (1,1) 
or 5*^ = (2,0) and S2 = (0,2). The first choice corresponds to a symmetric 
traceless tensor hfj,iy, while the second and the third choice are the irreducible 
pieces of the tensor R^v^pa-, which has the symmetry properties of the Reimann 
curvature tensor. Clearly, we will then take the (1, 1) irrep for the graviton, 
and one would then use the -ft'(i,i),2 propagator to construct the interaction 
vertices. 

The realistic interactions also require an introduction of the internal sym- 
metry group G. In the framework of the spin foam models, the simplest way 
to do this is to replace the category of representations Cat{G) with the cat- 
egory Cat{G X G). The matter irreps would be then given by {S, A5), where 
As is the corresponding irrep of G. The gravitational sector irreps would be 
(A^, 1), where 1 is the trivial irrep of G. The matter action would again have 
the form ( pSf ) in the topological gravity case, or (Q) in the non-topological 
gravity case, but now the products of the matter fields have to be contracted 
by the intertwiners C'^^'"^*' for the internal group irreps. 



8 Conclusions 

The constructions we have presented are inspired by the formalism of particle 
quantum field theories, and it is a straightforward extension of the existing 
gravitational spin-foam models. Given that the purely gravitational models 
can be related to a discretized path-integral quantization of the BF theory 
it is reasonable to expect that our model could also be related to an 



analogous discretized path-integral quantization of the BF theory with mat- 
ter. One would then need to rewrite the Einstein-Hilbert action with matter 
in an appropriate form, for example 



Si = J{BAR) + J {B AeA^-f{d + iu^)^) 

+ J{BAB)(^^r'U + g^rMa + ---) , (62) 

where u) is the spin connection, R = du + {uj A u) is the corresponding 
curvature two-form, Aa = e^A^, fab = e'^e'^f^u and the dots stand for other 
matter fields. (, ) is an appropriate contraction of the 5*0(3, 1) indices a, b, 
7 are the appropriate gamma matrices, /x, z/, ... are the curved space indices 
and are the inverse fierbeins. 
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Note that we have written the matter couphngs in terms of the tangent 
space components Aa, ... . These fields carry only the 5*0(3, 1) representation 
indices, and they would be the prime candidates for the continuous space 
hmit of our spin-network fields. This would also mean that one needs to 
introduce the fierbein form = e°;dx^ in the BF formalism for matter, since 
it is not possible to eliminate it from the action (^) by using the constraint 
B = e A e. However, in c? = 3 case there is no such problem, since B = e. 

As far as the graviton-like field h{x) is concerned, its continuous space 
analog would be hab{x) = e'^elh^u{x). In order to interpret this field as the 
graviton, one would need to implement the constraint h^y = g^y — gf^^, where 
gfiu = e^Ci^a is the spacetime metric, and g'^^ is a given background spacetime 
metric. This would then imply that the amplitude for a transition from 
a spatial spin net with n gravitons to a spatial spin net with m gravitons 
should not involve a sum over the background irreps N, since the background 
is fixed. That amplitude would be then given by the extension of the formula 
(18). Hence it is important to study further the model, in order to better 
understand the continuous spacetime limit, since there is a possibility that 
one would recover the perturbative gravity scattering amplitudes in a fixed 
background spacetime geometry. This also applies to the matter fields, so 
that this type of study will determine the correctness of our model. 

The fact that we have taken the matter spin net fields in the non-topological 
gravity model to have indices from the representations of the 5*0 (3) and 
SO (2) subgroups, means that the transition amplitudes will be 5*0 (3) or 
50(2) group covariant objects. One then wanders what happens with the 
Lorentz covariance. There is no such problem for the spin-half fermions, 
or for the (j, 0) and the (0,j) matter irreps. However, if we want to have 
the most general interactions, we need to include (j, k) irreps with jk ^ 0. 
Hence we encounter a spin-foam analogue of the particle field theory spinor- 
vector interaction problem. The only possible answer is that there should 
be a spin-foam version of the gauge invariance which reduces the Lorentz 
covariant expressions to 50(3) or 5*0(2) covariant gauge-fixed expressions. 
In that case the expressions ( |5^ ) and ( p5D can be understood as gauge-fixed 
actions. One can think of these actions as a spin-foam generalization of 
the light-cone gauge field theory actions |BD[. Actually, this analogy with 



the light-cone gauge field theory is even closer, since the finite-dimensional 
SO{n) irreps are unitary, and hence the corresponding spin-foam model in- 
volves only the unitary representations, in contrast to the topological gravity 
spin foam model, where the matter irreps are pseudo-unitary. 
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Obtaining the Lorentz covariant expressions in the non-topological grav- 
ity case would then amount to couphng simple unitary irreps to pseudo- 
unitary irreps S, and then forming the actions of the type (|53|). These actions 
should be invariant under a spin-foam version of the gauge transformations, 
so that one would obtain the action ( ^4|) in a particular gauge. A better 
understanding of these issues deserves a further study. 

Note that in the case of spin foam models where A2{A) is not dim A, it 
is not known what is the group field theory formulation. However, the field 
theory formulation is not necessary for the determination of the amplitudes. 
We used it in order to get a physical insight for our construction. Therefore, 
when A2 7^ dim A, we will postulate the amplitude formulas (|T7|) and (18), 
with Ae and A^ given by the same expressions as in the Af = dim A case. 

An important issue which has to be explored is the question of the per- 
turbative finiteness of the model. As in the purely gravitational case 



this will boil down to the study of the asymptotic behaviour of the spin-net 
propagators K]\r and Ks^s- Note that in the Euclidian case TrKs^s is the 
same as K]\r, where N = (s, s), so that one expects that the Euclidian spin- 
foam matter model should be finite. In the Lorentzian case there is no such 
a simple relation between the propagators, so that one needs to make a more 
detailed analysis. Alternatively, one can formulate a quantum group version 
of the construction we gave, which could be made automatically finite if the 
gravitational irreps are from a finite set, and if there are finitely many matter 
irreps. 

The issue of the semiclassical/continuous space limit of the spin foam 
models needs a more detailed study. There are no rigorous statements about 
the semiclassical theory. In the purely gravitational case one can show that in 
the limit of large angular momenta one obtains amplitudes which are sums of 
terms proportional to exp(±zS'/j), where Sr is the Regge form of the Einstein- 
Hilbert action Actually, because of the form of the propagator K^, see 
( pop and (^71), this seems to be a good approximation even for arbitrary 



values of the angular momenta. Hence there are good indications that one 
is on the right track. In the case of matter, one would then need to obtain 
expressions containing exp{±iSR ± iSgm), where Sgm is the simplical version 
of the matter actions This type of analysis could also settle the question 
of which fields are really massive and which ones are massless. 
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Figure I: A spin foam Fcynman diagram describing propagation of a matter 
field, which is represented by a dotted hne. 
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Figure II: The simplest spin net graph which gives a 7-vaIent vertex for a 
spin foam Feynman diagram. 




Figure III: A nontrivial spin net graph which gives a matter propagator for 
a spin foam Feynman diagram. 
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Figure IV: An example of a more complex spin net graph determining the 
7-valent spin foam vertex. 




Figure V: A more complex spin net graph determining the matter spin foam 
propagator. 
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Figure VI: Diagrammatic representation of the equation 
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